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1. INTRODUCTION 

Since its invention, NMR spectroscopy has developed 
from a technique for studying physical phenomena such 
as magnetism into a tool for acquiring information about 
molecules in chemistry and biology. Furthermore, it was 
pointed out early on (1955), almost as an anecdote, that 
nuclear spins could also be used for storing information 
using spin echoes 111 

This insight beautifully illustrated a notion that was 
developed in a very different context: information is 
physical and cannot exist without a physical represen- 
tation B In recent decades, the relationship between 
physics and information has been revisited from a new 
perspective: could the laws of physics play a role in how 
information is processed? The answer appears to be yes. 
If information is represented by systems governed by the 
laws of quantum mechanics, such as nuclear spins, an en- 
tirely new way of doing computation, quantum computa- 
tion (QC), becomes possible. Quantum computing is not 
just different or new; it offers an extraordinary promise, 
the capability of solving certain problems which are be- 
yond the reach of any machine relying on the classical 
laws of physics! 

The practical realization of quantum computers is still 
in its infancy. Interestingly, over 40 years after the initial 
suggestion of using spins to represent (classical) informa- 
tion, NMR spectroscopy has become the workhorse for 
experimental explorations of quantum information pro- 
cessing. 

In this article, we first explain how quantum computers 
work and why they could solve certain problems so much 
faster than any classical computer. Next, we describe 
how quantum computers can be implemented using NMR 
techniques and what is involved in designing and imple- 
menting QC pulse sequences, preparing a suitable initial 
state and interpreting the output spectra. We close with 
an overview of the state-of-the-art and the prospects for 
NMRQC and other QC implementations. 

Good reviews of quantum computation and informa- 
tion can be found in Refs. |jj-p| A good intro to and a 
specialized review of NMRQC are given in Ref. ^ and 
Ref. respectively. 



2. QUANTUM COMPUTATION 

2.1. INEPT as a computation 

We start from a familiar place for many NMR spectro- 
scopists, the INEPT pulse sequence (see INEPT). This 
sequence was designed to transfer polarization from a 
high 7 nucleus to a low 7 nucleus. However, it can also 
be viewed as a logic gate (Fig. which flips one spin 
conditioned upon the orientation of a neighbouring spin. 



If we arbitrarily assign "0" to a spin up and "1" to a 
spin down, we can think of spin-1/2 nuclei as bits in a 
digital computer. We remind the reader that bits ("0" 
or "1") can be used to represent numbers in a similar 
way to decimal numbers ( "0" through "9" ) . In decimal 
representation, 503 means 3 x 10° + x 10 1 + 5 x 10 2 ; 
similarly, in binary representation, 1101 means 1x2° + 
0x2 1 + lx2 2 + lx2 3 (which corresponds to the number 
13 in decimal representation). All computers today pro- 
cess information in binary representation, and are able 
to rapidly perform multiplication and division (repeated 
addition and substraction) as well as the most complex 
computations, via a sequence of very simple, elementary 
operations, called logic gates, acting on just one or two 
bits at a time. 

In this binary framework, the INEPT sequence corre- 
sponds to an elementary two-bit operation known as the 
controlled-NOT or for short CN OT g ate (within the phase 
corrections discussed in Section 3T ) . The CNOT gate per- 
forms a not operation on one bit, flipping it from "0" 
to "1" or from "1" to "0", if and only if the value of a 
second bit is "1" . The input to the logic gate is the ini- 
tial state of the spins, and the output is the final state of 
the spins. The four possible input values and the corre- 
sponding output values are tabulated in Fig. |^. 

The CNOT combined with single-spin rotations pro- 
vides for a universal set of logic gates. This means that 
any computational task can be implemented using a suf- 
ficiently large number of nuclear spins simply by con- 
catenating cnots and single-spin rotations in the proper 
way ETlJ. In summary, 

spin-1/2 nuclei in a molecule can serve as bits in a 
computer, and pulses and delay times provide universal 
logic gates. 



2.2. Quantum parallellism 

Everything we discussed so far was purely classical. 
The Bloch-sphere picture of Fig. |l| reinforces this classi- 
cal view of the spins; however, nuclear spins are really 
quantum objects. In Dirac notation, the state of a spin 
can be denoted |0) for a spin in the ground state (along 
z), and |1) for a spin in the excited state (along — z), 
corresponding to the two classical values for a bit ("0" 
and "1"). Now, a spin said to be 'along the x axis' is in 
reality in a coherent superposition state of spin up and 
spin down, written as (|0) + \l))/^/2, a spin 'along the y 
axis' is in the state (|0) + i|l))/-\/2, etcetera. A spin-1/2 
particle is thus more than just an ordinary bit. 

Any two-level quantum system, such as a spin-1/2 par- 
ticle, can serve as a quantum bit (qubit). 
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The difference between the quantum and classical de- 
scriptions becomes clear as soon as more than one quan- 
tum particle is considered. For example, it is well-known 
that the state of n interacting spins- 1/2 cannot be de- 
scribed simply by n sets of coordinates on the Bloch 
sphere. In order to include phenomena such as multiple- 
quantum coherence, we need recourse to 4™ — 1 real num- 
bers in the product operator expansion or equivalently to 
density matrices of dimension 2™ x 2™. Furthermore, the 
evolution of a closed system of n spins can only be de- 
scribed by 2™ x 2™ unitary matrices (see Liouville Equa- 
tion of Motion). The number of degrees of freedom 
that need to be specified in a classical description of the 
state and dynamics of n coupled spins thus increases ex- 
ponentially with the number of spins. 

Richard Feynman proposed in 1982 that the exponen- 
tial complexity of quantum systems might be put to good 
use to. simulate the dynamics of another quantum sys- 
tem lii, a task which requires exponential effort on a clas- 
sical computer. David Deutsch extended and formalized 
this idea in 19S5, and introduced the notion of "quantum 
parallellism" E3. 

Consider a (classical) logic gate which implements a 
function / with one input bit x and one output bit f(x). 
If x — 0, the gate will output /(0); if x = 1, the output 
will be /(l). The analogous quantum logic gate is de- 
scribed by a unitary operation which transforms a qubit 
as 

|0>~|/(0)) and |1>~|/(1)>. 

However, due to the possibility of preparing coherent su- 
perposition states and to the linearity of quantum me- 
chanics, the same gate also performs the transformation 



|0) + |1) 
V2 



1/(0)) + 
V2 



In this sense, it is possible to evaluate f(x) for both in- 
put values in one step! Next consider a different logic 
gate which implements a function g(x) with two input 
bits. We can prepare each qubit in a superposition of 
"0" and "1" . Formally, the state of the joint system is 
then written as (|0) + |1)) <g> (|0) + |l>)/2. Leaving out 
the tensor product symbol as well as any normalization 
factors, the state can be written as (|0) + |1))(|0) + |1)), or 
|0) |0) + 10) 1 1) + 1 1) |0) + 1 1) 1 1) , which is further abbreviated 
to 1 00) + |01) + 1 10) + |11). Therefore, a set of two spins 
can be in a superposition of the four states "00" , "01" , 
"10" and "11". A quantum logic gate implementing g(x) 
then transforms this state as 

|00)+|01)+|10)+|11) » | ff (00))+|s(01))+|ff(10))+|«/(ll)) 

Thus the function has been evaluated for the four pos- 
sible input values in parallel. In general, a function of 
n qubits implemented on a quantum computer can be 



evaluated for all 2™ input values in parallel! In contrast 
to classical computers, for which the number of parallel 
function evaluations increases at best linearly with their 
size, 

the number of parallel function evaluations grows ex- 
ponentially with the size of the quantum computer (the 
number of qubits). 

Obviously, this is true only so long as the coherent su- 
perposition states are preserved throughout the compu- 
tation. This means that the computation should be com- 
pleted before quantum coherence is lost due to "decoher- 
ence" processes (in NMR spin-spin and spin-lattice relax- 
ation; see Relaxation: An Introduction). Since some 
degree of decoherence is unavoidable, practical quantum 
computers appeared virtually impossible to build, until 
quantum error correction was conceived, as discussed in 



Section 2.4 



Even if the coherence time is long compared to the 
duration of a typical logic gate and quantum error cor- 
rection is employed, can we really access the exponential 
power exhibited by quantum systems? The postulates of 
quantum mechanics dictate that a measurement — which 
induces instantaneous and complete decoherence — of a 
qubit in a superposition state |/(0)) + |/(1)) will give ei- 
ther "f(0)" or "f(l)", with equal probabilities. Similarly, 
after doing 2™ computations all at once, resulting in a 
superposition of 2 n output values, measurement of the 
quantum bits randomly returns a single output value. A 
more clever approach is thus needed: 

exploiting quantum parallellism requires the use of 
quantum algorithms. 



2.3. Quantum algorithms 

Special quantum algorithms allow one to take advan- 
tage of quantum parallellism in order to solve certain 
problems in far fewer steps than is possible classically. 
When comparing the capability of two computers to solve 
a certain type of problem, 

the relevant criterion is not so much what resources 
(time, size, signal-to-noise ratio, . . .) are required to solve 
a specific instance of the problem but rather how quickly 
the required resources grow with the problem size. 

A particulary important criterion is whether the re- 
quired resources increase exponentially or polynomially 
with the problem size. Exponentially difficult problems 
are considered intractable — they become simply impos- 
sible to solve when the problem size is large. In contrast, 
polynomially difficult problems are considered tractable 
or possible to solve. The interest in quantum computing 
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is based on the fact that 

certain problems which appear intractable (resources 
grow exponentially with problem size) on any classical 
computer are tractable on a quantum computer. 

This was shown in 1994 by Peter Shor, almost 10 years 
after Deutsch introduced quantum parallcllism. Shor's 
quantum algorithm £3 allows one to find the period of 
a function exponentially faster than any classical algo- 
rithm. The importance of period-finding lies in that it 
can be translated, using some results from number the- 
ory, to finding the prime factors of integer numbers, and 
thus also to breaking widely used cryptographic codes. 
These codes are based precisely on the fact that no effi- 
cient classical algorithm is known for period-finding or 
factoring, i.e. the effort required to factor a number 
on classical computers increases exponentially with the 
number of digits of the integer to be factored. In contrast, 
Shor's algorithm is efficient: on a quantum computer, the 
effort to factor an integer increases only polynomially (to 
be precise, the third power) with the number of digits 
of the integer. As a result, while factoring a 1000-digit 
number is believed to be beyond the reach of any ma- 
chine relying on the classical laws of physics, such a feat 
could be accomplished on a quantum computer. 

The first^niantum algorithm was invented by Deutsch 
and Jozsa tfl (1992). This algorithm allows a quantum 
computer to solve with certainty an artificial mathemat- 
ical problem known as Deutsch's problem. Even though 
this algorithm does not have much application, it is his- 
torically significant as it provided the first steps towards 
Shor's algorithm, and because it is a simple quantum al- 
gorithm that can be experimentally tested. 

Another class of quantum algorithms-was discovered in 
1996 by Lov Grover. These algorithms ll^l allow quadratic 
speed-ups of certain search problems, for which there is 
no better approach classically than to try all L candi- 
date solutions one at a time. A quantum computer using 
Grover's algorithm needs to make only \J~L such trials. 
Even though this speed-up is only quadratic rather than 
exponential, it is still significant. 

The last currently known application of quantum_cpm- 
puters lies in simulating other quantum systems E3, as 
Feynman conjectured. Even a computer consisting of 
no more than a few dozen quantum bits could outper- 
form the fastest classical computers in solving important 
physics pxablems, such as calculating the energy levels of 
an atom O. 

In the remainder of this section, we will briefly review 
the structure of Shor's algorithm, because it is so impor- 
tant and at the same time gives good insight into how 
quantum computing works (for a more detailed explana- 
tion, see Refs. |^ and [l3]). The crucial step in Shor's 
factoring algorithm is the use of the quantum Fourier 



transform (QFT) to find the period r of the function 
f(x) = a I modM, which means f(x) is the remainder af- 
ter division of a x by M, where M is the integer to be 
factored, and fLJS. an integer which is more or less ran- 
domly chosen oE3. 

The QFT performs the same transformation as the 
(classical) fast Fourier transform (FFT), but can be com- 
puted exponentially faster. As always, we don't have ac- 
cess to all the individual output values; the QFT merely 
allows us to sample the FFT but as we will see this suffices 
for period-finding. The FFTjv takes as input a string of 
N complex numbers Xj and produces as output another 
string of N complex numbers yk, such that 

w = ^£V**'". (i) 

For an input string with numbers which repeat them- 
selves with period r, the FFTat produces an output string 
with period N/r, as illustrated in the following four ex- 
amples for N = 8 (the output strings have been rescaled 
for clarity) 

r input string output string N/r 

8 1 h 11111111 1 (a) 
4 1 1 ^ 10101010 2 (b) 
2 10101010 h 1 1 4 (c) 
1 11111111 w 1 8 (d) 

If r does not divide N, the inversion of the period is 
approximate. Furthermore, the FFT turns shifts in the 
locations of the numbers in the input string into phase 
factors in front of the numbers in the output string: 

10001000 h-> 10101010 (c) 

1 1 i — ► 1 -i -1 i (f) 

1 1 h-> 10-1010-10 (g) 

1 1 i — * 1 i -1 -i (h) 

The QFT performs exactly the same transformation, 
but differs from the FFT in that the complex numbers 
are stored in the amplitude and phase of the terms in 
a superposition state. The amplitude of the 1 000) term 
represents the first complex number, the amplitude of 
the |001) term the second number and so forth. For 
clarity, we will label the states |000), |001), . . . | 111) as 
|0), |1), . . . |7). As an example, we see from (f) that 
the QFT transforms the state |1) + |5) to the state 
|0> — i|2> — |4> 

The QFT is incorporated in actual quantum algo- 
rithms as outlined in Fig. |^. A first register (group of 
qubits) is prepared in a superposition of all its possible 
states. A second register is initialized to the ground state 
(for factoring a number M, the size of the second regis- 
ter must be at least log 2 M and the first register must 
be at least twice as large). For example, if register 1 has 
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three qubits and register 2 has two qubits, the state of 
the system is prepared in 

(|0) + |1) + |2) + |3) + |4) + |5) + |6) + |7))|0). (2) 



Then the function f(x) is evaluated (in NMR by apply- 
ing a pulse sequence, as discussed in Sections 2.1 and 3.1 ) 



where x is the value of the first register, and the output 
value f(x) is stored in the second register. Since the first 
register is in an equal superposition of all \x) , the func- 
tion is evaluated for all values of x from to 7 in parallel. 
For example, let f(x) — 3 for even x and f(x) = 1 for 
odd x, which means the period r is 2 (in real applica- 
tions, we have a description of / but do not know r in 
advance). Evaluation of f(x) then transforms the state 
of eq. H to the state 

|0)|3) + |1}|1} + |2)|3) + |3)|1) 

+ |4)|3) + |5>|1) + |6)|3) + |7>|1> (3) 

= (|0) + |2) + |4> + |6))|3) + (|1) + |3) + |5) + |7»|1) . 

(4) 

We pause to point out that this state is entangled, which 
means that it cannot be written as a product of single- 
qubit states. The state |00) + |01) + 1 10) + |11) is an ex- 
ample of an unentangled state, because it can be written 
as (|0) + |1))(|0) + |1)), a product of single-qubit states. 
In contrast, 1 00) + |11) is a simple example of an entan- 
gled state. Entanglement has no classical analogue and 
is believed to lie at the heart of the exponential speed-up 
offered by quantum computation. 

In order to appreciate the role of the QFT, suppose we 
now measure the second register in Eq. || (this measure- 
ment can be left out but simplifies the explanation). The 
state of the first register will collapse to either 

|0) + |2) + |4) + |6) or |1) + | 3 ) + |5) + |7), (5) 

depending on whether the measurement of register 2 gave 
"3" or "1". We see that all eight possible outcomes are 
still equally likely, so the result of a measurement does 
not give us any useful information. But if we apply the 
QFT, the first register will be transformed to 



|0) + |4) or |0) - |4) 



(6) 



Now a measurement of the first register does give useful 
information, because only multiples of N/r are possible 
outcomes, in this example "0" and "4" . 

This concludes the quantum part of the computation. 
From the measurement result, a classical computer can 
efficiently calculate the inverted period N/r, and thus 
also r, with high probability of success using results from 
number theory. Now that r is known, the factors of the 
integer M can be quickly computed as well, with high 



probability (the probability of success can be further in- 
creased by repeating the whole procedure a few times). 

We close with two final remarks on quantum algo- 
rithms. (1) Quantum computing cannot offer any speed- 
up for many common tasks, such as adding up two num- 
bers or word processing, which can already be completed 
efficiently on a classical computer. (2) There are many 
exponentially difficult problems which no currently avail- 
able quantum algorithm could help solve faster than is 
possible classically. It would be somewhat disappoint- 
ing from a practical viewpoint if no other applications 
were found; however, our understanding of the connec- 
tion between physics and information and computation 
has already changed dramatically. 



2.4. Quantum error correction 

Any quantum computation must be completed within 
the cohere nce time, in NMR T2 and T\, as pointed out 



in Section 2.2. T\ and T2 processes alter the state of the 
qubits and are therefore a source of errors. For many 
years, this requirement led to wide-spread pessimism 
about the practicality of quantum computers. In 1995, 
however, Peter Shor and Andrew Steane_independently 
discovered quantum error correction l 46 ^ 7 ! and showed 
that 

It is possible to correct for truly random errors caused 
by decoherence. 

This came as quite a surprise, because quantum error 
correction had to overcome three important obstacles: 
(1) the no-cloning theorem, which states .that it is not 
possible to copy unknown quantum states 13, (2) measur- 
ing a quantum system affects its state and (3) errors on 
qubits can be arbitrary rotations in Hilbert space, com- 
pared with simple bit flips in classical computers. Quan- 
tum error correction requires many extra operations and 
extra qubits (ancillae), though, which might introduce 
more errors than are corrected, especially because the ef- 
fect of decoherence increases exponentially with the num- 
ber of entangled qubits, much in the same manner that 
multiple quantum coherences decay exponentially faster 
than single quantum coherences. Therefore, a second sur- 
prising result cB was that 

provided the error rate (probability of error per elemen- 
tary operation) is below a certain threshold, and given a 
fresh supply of ancilla qubits in the ground state, it is pos- 
sible to perform arbitrarily long quantum computations. 

The threshold error rate is currently estimated Efl to 
be about 0.001%. The actual error rate in NMRQC is ap- 
proximated by 1/2 JT2, where 2JTi is roughly the num- 
ber of operations that can be computed within the coher- 
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ence time. For small molecules the error rate is typically 
on the order of 0.1% to 1%, two to three orders of magni- 
tude too high. But a remarkable implication of quantum 
error correction is that if (1) a molecule is found which 
achieves the accuracy threshold and (2) select spins can 
be fully polarized, both T\ and T2 could in principle be in- 
finitely lenghtened by applying an error correction pulse 
sequence. 



3. NMR QUANTUM COMPUTERS 

3.1. Pulse sequence design 

The translation of abstract quantum algorithms or 
function evaluations into actual pulse sequences may ap- 
pear obscure at first sight. However, systematic tech- 
niques exist to make pulse sequence design rela- 
tively straightforward. The starting point is that 

each quantum algorithm can be described by a sequence 
of transformations under unitary operators. 

Such unitary transformations represent rotations in 
Hilbert space (a multidimensional extension of the Bloch- 
sphere). Examples of unitary transformations are evolu- 
tion during RF pulses and free evolution under the sys- 
tem Hamiltonian; relaxation processes give rise to non- 
unitary transformations. Once the desired unitary oper- 
ators have been identified, 

arbitrary unitary operators can be translated into se- 
quences of single-qubit rotations and CNOT gates. 

These building blocks can be readily implemented in 
Decompositions into other sets 
}d can be helpful 
In any case, it 



NMR (see Section p.2| ). 
of elementary gates are also possibL 
in simplifying the pulse sequences 
is crucial that the duration of the pulse sequence design 
process as well as the length of the resulting pulse se- 
quence not increase exponentially with the problem size. 

We now point out two important distinctions between 
QC and conventional pulse sequences. On the one hand, 
QC sequences must be more general: 

QC sequences must perform the desired transformation 
for arbitrary input states. 

In contrast, conventional sequences are often designed 
assuming a particular input state. As a first example of 
this difference, the sequence of Fig. [TJ assumes that both 
spins are in Zeeman states, i.e. aligned along ±z. It 
implements the unitary operator 



INEPT 



/l CP 

i 

1 

V0 -i 0, 



(7) 



which is similar to but different from the unitary operator 
for the CNOT gate, defined as 



C^CNOT — 




(8) 



implemented, for example, by 90° 90^ z 90^ 1/2 J ab 90 b _ y - 
As a second example, consider the so-called Hadamard 
gate, defined as 



Had 



Wl 1 

V2\l -1 



(9) 



This gate creates a superposition state starting from a 
basis state: it transforms |0) to |0) + |1) (z to x in the 
Bloch sphere) and |1) to |0) — |1) {—z to —x). At first 
sight, this transformation could be done simply via a 90 a 
pulse. However, the unitary operator for 90j, 



U$a.. — 



= ±( 1 



1 



V2 V 1 1 



(10) 



is different from C% a( j; e.g. applying t^TTad ^ w i ce nas 
no net effect, but applying [/90 twice produces Uiso ■ A 
possible sequence to implement £% ac j exactly is 90yl80 x . 

On the other hand, QC sequences can be more specific: 

QC sequences can be specialized for a specific molecule 
using full knowledge of its spectral properties. 

In contrast, conventional sequences must work for any 
molecule, because the spectral properties of the molecule 
are usually not known in advance. Exact knowledge of 
the chemical shifts and J-coupling constants allows one 
not only to greatly simplify the pulse sequences, but also 
to achieve much more accurate unitary transformations 
than would otherwise be possible. 

Finally, while systematic procedures exist to design a 
pulse sequence, there is a need to develop tools to find 
the pulse sequence of the shortest duration and with the 
fewest RF pulses. Even small-scale quantum computa- 
tions easily involve tens to hundreds of gates acting on 
multiple spins and precise control of the spin dynamics is 
difficult to maintain throughout such long sequences of 
operations, as shown in the next section. 



3.2. Implementation of computations 

The implementation of quantum computations with 
NMR can be based on single-spin rotations and CNOT 
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gates, since any quantum algorithm can be translated 
into these building blocks. Although these elementary 
operations appear quite easy to implement, 

the requirements for precision in QC experiments are 
unusually high, due to the large number of pulses and the 
quantitive nature of the information contained in the out- 
put spectra. 

Implementation of accurate single-spin rotations about 
an axis in the x-y plane is relatively easy in heteronuclear 
molecules; yet, it can be very demanding for homonu- 
clear spin systems because spin selectivity requires longer 
pulses resulting in considerable coupled evolution of the 
spins during the pulses Eij. Clearly, some degree of 
homonuclearity is unavoidable when more than a handful 
of qubits are involved. 

We therefore begin by reviewing the requirements 
for pulse shaping E3 (see Shaped Pulses; Selective 
Pulses). First, the magnetization corresponding to each 
of the lines in a multiplet must be rotated about ex- 
actly the same axis and over exactly the same angle, 
i.e. off-resonance effects due to line splitting must be re- 
moved. This requires self-refocusing shaped pulses or tai- 
lored composite pulses E3 (see Composite Pulses). Sec- 
ond, the effect of J-couplings between unselected spins 
must be removed, either during the pulse or later in the 
pulse sequence. Third, all pulses must be universal ro- 
tors, i.e. the rotation must be independent of the initial 
state of the spin. Fourth, the unselected spins must not 
be affected by the RF irradiation. This last requirement 
is difficult to satisfy because of transient Bloch-Siegert 
effects £3 which can result in substantial (tens of de- 
grees) phase shifts of nearby unselected spins. However, 
it is possible-tOpestimate and compensate for the Bloch- 
Siegert shift E3'E3. Finally, simultaneous (as opposed to 
consecutive) pulses at two or more nearby frequencies 
are desirable in order to keep pulse sequences short, but 
transient Bloch-SiegepLsMfts greatly deteriorate such si- 
multaneous rotations EZpEM Nevertheless, accurate simul- 
taneous rotations at nearby frequencies can be achieved 
using a special correction technique E£l Still, simulta- 
neous pulses on welLcoupled spins may excite multiple 
quantum coherences l52l 

There are a number of hardware requirements for suc- 
cessful execution of QC experiments. Good RF coil ho- 
mogeneity is crucial in avoiding excessive signal attenua- 
tion and related errors. Furthermore, it is desirable that 
one frequency source and transmitter board be available 
per qubit. If there arc more qubits than spectrometer 
channels, the carrier frequency must be jumped to the 
appropriate frequencies throughout the pulse sequence, 
or phase ramping techniques must be employed c3. A 
dedicated frequency source for each qubit also makes it 
easy to keep track of the rotating frame of each spin and 
apply all the pulses on any given spin with the correct 



relative phase. Xhis removes the need to refocus chemical 
shift evolution E3, which involves extra pulses. Alterna- 
tively, software rotating frames can be created by detailed 
bookkeeping of the time elapsed since the beginning of 
the pulse sequence such that the evolution of the rotating 
frame of any given spin with respect to the carrier refer- 
ence frame can be calculated. The phases of the pulses 
throughout the pulse sequence, as well as t he receiver 
phase, can then be adjusted accordingly P 6 M 2 4 The same 
technique can be used to implicitly realize single-spin ro- 
tations about z. Alternatively, z-rotations can be imple- 
mented using resonance offsets or composite pulses E 2 ! 

Two strategies exist for implementing CNOT gates 
(both assume first order spectra). If all the spins are mu- 
tually coupled, CNOT's can be realized via line-selectisze 
pulses which invert specific lines within a multiplet 113. 
In practice, it is usually more convxmiipt to use pulse 
sequences such as the one in Fig. |l| For molecules 

with several coupled spins, the sequence of Fig. [j] must 
be expanded with extra pulses to refocus the undesired 
J-couplings; systematic methods exist to design good re- 
focusing schemes E2n2il'E3. A CNOT between two uncoiL- 
pled spins can be realized by swapping qubit states MJ 35 ! 
For example, for a CNOT between two spins a and c which 
are not mutually coupled but which are both coupled to 
a third spin 6, the procedure is as follows: apply a pulse 
sequence which swaps the state of a and b (via CNOT {, 
CNOTb a CNOT a b), then perform a CNOT between b and c, 
and then swap a and b again. The net result is CNOT ac ; 
spin b is unaffected. It is thus not necessary that all spins 
be pairwise coupled as long as the network of couplings 
includes all n spins. 

An alternative to imposing the correct evolution on 
all the spins at all times, both for RF pulses and for 
CNOT-type gates, is to allow erroneous evolutions which 
will later be reversed. Such techniques have been highly 
successful in certain standard pulse sequences (see De- 
coupling Methods), but are much harder to develop 
for the non-intuitive and non-transparent QC pulse se- 
quences. Nevertheless, it has been shown experimentally 
that a large degree of cancellation of erroneous evolutions 
is possible even in QC experiments: about 300 two-qubit 
gates involving over_1350 RF pulses have been success- 
fully concatenated E3. A general methodology designed 
to take advantage of this possibility has yet to be devel- 
oped. 

From this discussion, it will be clear that 

the selection of a suitable molecule is crucial for 
NMRQC. 

The desired properties are (1) sufficiently large chemi- 
cal shifts for good addressability, (2) large coupling con- 
stants, while maintaining first-order spectra, for fast two- 
qubit gates (or a coupling network that matches the pat- 
tern of connectivities needed for the algorithm), and (3) 
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long T2's and Ti's in order to allow time to execute many 
logic gates. Furthermore, high-7 nuclei are desirable for 
good sensitivity. More mundane but equally important 
requirements are that the molecule be stable, synthesiz- 
able, soluble and safe. 

3.3. State initialization 

Apart from experiments designed to produce non- 
thermal spin polarizations, setting up a proper initial 
state for the nuclear spins is a concept worth revisiting 
for NMR spectroscopists. Since this is a crucial step in 
quantum computing, this entire section is therefore de- 
voted to state initialization. 

Most quantum computations require a pure initial 
state, for example a set of fully polarized spins, in the 
state 1 00 . . . 0). However, nuclear spins in thermal equilib- 
rium at room temperature are in an almost fully random 
state: for typical magnetic field strengths, the ground 
(|0)) and excited (|1)) state probabilities differ by only 
about 1 part in 10 5 . The spins are then said to be in 
a mixed (non-pure) state. The polarization could be 
increased using hyperpolarization techniques (see Un- 
usual signal enhancement: optical pumping & hy- 
perpolarized inert gases and Optically enhanced 
magnetic resonance) but the state of the art is still 
very far from cooling nuclear spins into the ground state. 

The conceptional breakthrough which made NMR 
quantum computation possible at room temperature 
was the^the concept of effective pure or pseudo-pure 
states US 

effective pure states are mixed states which produce the 
same signal as a pure state to within a scaling factor. 

The signature of an effective pure state for n spins is 
that all but one of the 2™ populations are equal, and 
that no coherences arc present. The density matrix then 
consists of an identity component and a pure state com- 
ponent. The identity density matrix is not observable in 
NMR since only population differences are observed, and 
furthermore does not transform under unitary evolutions 
(UIU' S — I). Thus the visible signal is produced solely 
by the one distinct population, corresponding to a pure 
state. In product operator notation EJ, the effective pure 
ground state is proportional to I x + S z + 2I Z S Z for two 
spins, to I z + S z + R z + 2I Z S Z + 2I Z R Z + 2S Z R Z + AI Z S Z R Z 
for three spins, and so forth. Figure |] shows how the ef- 
fective pure state preparation is manifest in the spectrum 
of one of 5 coupled spins. Characteristic of effective pure 
(basis) states is that only one line survives in each mul- 
tiplet. 

Three methods are known for preparing effective pure 
states starting from thermal equilibrium. 



(1) Logical labeling consists of applying a pulse 
sequence which rearranges the thermal populations such 
that a subset of the spins is in an effective pure state, con- 
ditioned upon the state of the remaining spins. Then the 
computation is carried out within this embedded sub- 
system e3. For example, the Boltzman populations for 
the states {|000), |001), |010>, |011), |100), 1 101) , |110), 
I 111)} for a homonuclear three-spin system deviate from 
the uniform background by {3a, a, a, —a, a, —a, —a, —3a} 
respectively, where a = 2 ^ T <C 1. After rear- 
ranging the populations for the eight spin states as 
{3a, —a, —a, —a, a, a, a, —3a}, the last two qubits are in 
an effective pure state conditioned upon the first qubit 
being |0) . As the total number of qubits n in the molecule 
increases, the relative fraction of effective pure qubits 
goes to 1, but the preparation sequence becomes complex 
quite rapidly for large n and the signal strength scales as 
n/2 n . ^ 

(2) Temporal averaging l!2l is similar to phase cy- 
cling (see Phase Cycling), since it consists of adding 
up the spectra of multiple experiments. However, in- 
stead of changing just the phase of some pulses, each 
experiment starts off with a different state prepara- 
tion sequence which permutes the populations. For 
two heteronuclear spins, adding together three exper- 
iments which yield respective population deviations 
{a, 6, —6, —a}, {a, —b, —a, 6} and {a, —a, 6, — 0} is equiva- 
lent to performing an experiment with population devia- 
tions {3a, —a, —a, —a}. For arbitary n, at least (2 n — l)/n 
experiments are needed E3, since the effective pure state 
is made up of 2™ — 1 product operator terms and the 
starting state, thermal equilibrium, contains n terms. 

(3) Spatial averaging li£| uses a pulse sequence contain- 
ing magnetic field gradients (see Field Gradients & 
Their Application) to equalize all the populations but 
the ground state population. Only one experiment is 
involved, but the preparation sequence quickly becomes 
unwieldy for large spin systems and the signal strength 
decreases exponentially with n. 

To date, temporal and spatial averaging have been the 
most popular choices for preparing effective pure states. 
Several hybrid schemes c3'E3 have also been developed 
which trade off complexity of the preparation steps for 
the number of experiments. Nonetheless, all these state 
preparation schemes have in common that 

creating effective pure states incurs an exponential cost 
either in the signal strength or in the number of experi- 
ments involved. 

Such an exponential overhead is of course not tolerable 
for quantum computations. The reason for this cost is 
that effective state preparation techniques simply select 
out the signal from the ground state population present 
in thermal equilibrium, and the fraction of the molecules 
in the ground state is proportional to n/2 n . computation. 
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A significant breakthrough by Schulman and Vazirani 
resulted in a method to cool a subset of the spins in 
a molecule down to the ground state without any ex- 
ponential overhead l 41 H 42 l The idea is to redistribute 
the entropy of the spins so that some have zero entropy 
(pure state) while the entropy of the remaining spins in- 
creases. In a sense, this method is an advanced polar- 
ization transfer scheme. Surprisingly both the length of 
the cooling pulse sequence and the number of spins that 
must be sacrificed increase only about linearly with the 
number of pure spins desired. Furthermore, the cooling 
algorithm approaches the entropic bound in the limit of 
large numbers of spins. However, with initial polarization 
a <C 1, a molecule with at least k/a 2 spins is needed to 
obtain k pure spins. This method is therefore impracti- 
cal when starting from equilibrium at room temperature, 
with a « 10~ 5 . Nevertheless, a combination of hyperpo- 
larization techniques and the Schulman- Vazirani scheme 
may some day become practical. In any case, despite the 
exponential cost incurred when preparing effective pure 
states, 

the highly random initial state represents no fundamen- 
tal obstacle to scalable quantum computation. 



3.4. Read-out 

Traditionally in NMR spectroscopy, only one nucleus 
is observed. In QC, however the concept of a single ob- 
serve channel and one or more decoupler channels does 
not apply: 

the output of a quantum computation is the final state 
of one or several spins. The final states of each of the 
output spins must thus be read out. 

If each of the output spins ends up in |0) or |1) (or 
in reality in the effective pure state corresponding to |0) 
or |1)), the answer can be read-out directly by applying 
a pulse which rotates the spin from ±z to ±x. With 
properly referenced receiver phase settings, the spectrum 
for each output spin then consists of either absorption or 
emission lines, indicating whether the output value of the 
corresponding bit is "0" or "1" (Fig. |||). 

If the output state is a superposition state, the situa- 
tion is a bit more complicated. For a single (as opposed 
to an ensemble) quantum computer subject to a "hard" 
measurement (assumed in Section 2.2), the superposition 
"collapses" to one of the terms in the superposition, with 
probabilities given by the square of the amplitude of each 
term. In contrast, 



measurements in NMR give a (bit-wise) ensemble av- 
eraged read-out. 



The output state of eq. || serves as an example: half 
of the molecules in the ensemble collapse to |0) ( 1 000) ) , 
while the other half collapses to |4) ( 1 100) ) . In other 
words, spins 2 and 3 always end up in "0" so their spectral 
lines are absorptive; in contrast, the signal of spin 1 aver- 
ages to zero because there are equally many molecules in 
which spin 1 ends up in "0" as in "1" . It is not clear that 
such bit-wise averages of probabilistic output states are 
generally sufficient to solve the problem of interest. For 
Shor's period-finding algorithm, this problem can be cir- 
cumvented □ by performing the classical post-processing 
steps (Section |2.3| ) on the quantum computer using some 
ancillae qubits — any classical-computation can also be 
done on a quantum computer In this way, the output 
state becomes the period r for all the molecules in the en- 
semble (as opposed to an average over all the multiples of 
N/r) and the measurement result becomes deterministic. 

Instead of recording the signal of each of the output 
spins, it is sometimes possible to use the extra informa- 
tion provided by the line splittings due to J couplings to 
derive the output state of several qubits from the spec- 
trum of a single spin. Since each of the lines in the mul- 
tiplet can be identified with specific states of the other 
spins (as in Fig. |J), the presence or absence of each line 
in the multiplet gives information about the state of the 
other spins (Fig. ||). 

Finally, while the spectra of a few select spins suffice 
to obtain the answer to a computation, the full density 
matrix conveys much more information. This extra in- 
formation can be used to expose the presence of errors 
such as multiple quantum coherences not visible in the 
single output spectra and furthermore is a useful tool 
for debugging pulse sequences. The procedure for re- 
constructing r the density matrix is called quantum state 
tomography c3c3. It consists of repeating the computa- 
tion multiple times, each time looking at the final state of 
the spins after applying different sets of read-out pulses 
which rotate different elements of the density matrix to 
observable positions. However, since this procedure in- 
volves on the order of 4™ experiments, it is practical only 
for experiments involving a few spins. 



3.5. State-of-the-art and outlook 

To date, only very simple demonstrations of quantum 
algorithms, all using liquid-state NMR techniques, have 
been carried out. Variations of C jroy cr's algorithm haw 
been demonstrated with two rM' l 49 !! 50 ! and three qubits E3, 
the Deulach-Jozsa algorithm with two c3'EiJ, three E3 
and five qubits, and the period-finding algorithm with 
five qubits A^niantum sinrulation has been imple- 
mented with two E3 and three E3 qubit s. a nd quantum 
error detection and correction with two and three I 57 ! 



qubits. Fina 
and observed 



a 7-spin coherence has been created 
(see Ref. |^ for additional references). 
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While these experiments indeed demonstrate the princi- 
ples of quantum information processing, they all involve 
far fewer qubits than would be needed to solve a problem 
beyond the reach of classical machines. 

Despite the rapid progress in recent years, scaling liq- 
uid state NMRQC to tens or hundreds of qubits may 
be impractical for several reasons, although none of 
them appear fundamental. In particular, as the num- 
ber of qubits increases: (1) the strength of the signal 
selected with current^tate initialization techniques de- 
creases exponentially Era'EjJ but Schulman-Vazirani cool- 
ing cil does not suffer any such exponential overhead (Sec- 
tion [j^ ); (2) the chemical shift separations unavoidably 
become smaller, but Lloyd showed that for universal 
quantum computation it suffices to have a linear chain 
d — abc — abc — ... — abc, in which only nearest neighbour 
spins are coupled and with only four distinct chemical 
shifts 8 a , 8b, 8 C and 84 ; (3) J-couplings become smaller or 
even unresolved, but this can 
ping spin states (Section |3~^ ) 
are not fundamental, the solutions all make the pulse 
sequences much longer. This would require increasingly 
longer T2's and Ti's for larger molecules, while in prac- 
tice the T2 , s and Ti's tend to become shorter (see Re- 
laxation: An Introduction). 

NMRQC has also brought up new theoretical issues. 
(1) Since only ensemble averaged results are available be- 
cause of the large number of molecules in a sample tube, 
some information is lost that would be available in an 
idealized quantum computer such as a single molecule 
at zero degrees Kelvin. For the known quantum algo- 
rithms, this information can be retrieved by performing 
classical post-processing steps on the quantum computer 
(Section 2.3). (2) Since the density matrix of nuclear 



circumvented by swap- 
32 l While these obstacles 



spins at room temperature is very close to the identity 
matrix, it is not possible to produce genuinely entangled 
states between the nuclear spins LtLsmall thermally polar- 
ized molecules in liquid solution E3. This observation has 
sparked a stimulating debate about the "quantumness" 
of NMR, because it implies that each of the states pro- 
duced in NMRQC experiments so far is classical. How- 
ever, all attempts to describe the dynamics of a set of 
coupled spins by an efficient classical model have been 
unsuccesful. It is thus conjectured that even though the 
states are classical, the dynamics of the spins is truly 
quantum mechanical E3, a proposition which will appear 
obvious to most NMR spectroscopists. In fact, it has 
also been the starting point of this introduction to NMR 
quantum computing. 



4. SUMMARY AND CONCLUSIONS 

In many respects, liquid state NMR provides an ideal 
test bed for elementary quantum computations. The de- 
gree of control over the evolution of multiple coupled 



qubits — the result of 50 years of technology develop- 
ment — the long relaxation times of nuclear spins and a 
set of new insights Era made it possible to perform certain 
computations in fewer steps than is possible using any 
classical machine. This is in itself a remarkable achieve- 
ment. 

It is unlikely that liquid state NMR could ever be used 
to solve problems faster than any classical machine, but 
it has already inspired many other NMR based propos- 
als for quantum computing. Liquid-crystal solvents have 
been used to partly reintroduce dipole-dipole couplings 
(see Liquid Crystals: General Considerations) to 
speed up the gate time and increase the number of gates 
possible within the coherence time til New molecular 
architectures based on liquid crystal solutions are now 
being investigated. Solid-state NMR near zero Kelvin 
could circumvent the state initialization problem, but 
also poses new questions in terms of addressability and 
coherence times (see Internal Spin Interactions & 
Rotations in Solids). Several approaches to solve these 
issues have been proposed E3'E3. Another proposal which 
received much attention consists of doing NMR on impu- 
rity atoms placed in a linear or two-dimensional array, 
with chemical shifts and couplings controlled by elec- 
trodes -placed on top of and in between the impurity 
atoms l£f . 

Furthermore, there is a plethora of very different ex- 
perimental approaches to building quantum computers 
(for an extensive review, see Ref . pa) . Four trapped 
ions have recently been entangled l££Tand a single quan- 
tum logic gate has been implemented between two pho- 
tons coupled with each other via interactions with an 
atom in a cavity. In the long run, the most scalable ap- 
proaches may be those based on solid-state technology, 
such as electron spins in quantum dots or magnetic fluxes 
in SQUIDs E3. 

It is clear that none of these proposals will be easy 
to implement since they all require substantial and in- 
novative development of technology. The success of any 
approach will depend on the ratio of the coherence time 
to the gate duration, i.e. how many gates can be com- 
pleted within the coherence time, and on the achievable 
degree of quantum control. 

Many of the problems in quantum control are simi- 
lar for different experimental systems, and we therefore 
expect that other quantum computer implementations 
will benefit from the ideas, concepts and solutions which 
arise from liquid state NMR experiments. In addition, we 
hope that some of the techniques developed within the 
context of quantum computation may find more general 
application in NMR. 

The possible payoff for successful quantum computing 
is tremendous: to solve problems beyond the reach of any 
classical computer. It is not clear at this point whether 
quantum computers will fulfill this promise, but in any 
case quantum computing has already provided an excit- 
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ing new perspective on NMR and, more broadly, on the 
connection between physics, information and computa- 
tion. 
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FIG. 1. The evolution of one of two coupled heteronuclear 
spins during an INEPT type pulse sequence, when the other 
spin is up (solid line) or down (dashed line). The rotating 
frame is set on resonance with the first spin so there is no 
need to refocus chemical shift. The usual read-out pulse is 
left out. The same pulse sequence can be applied to two 
homonuclear spins using spin-selective pulses. 



FIG. 5. Output spectra of the proton (left) and carbon 
(right) spins of 13 CHCl3 (disolved in a liquid crystal solution) 
for four different executions of Grover's search algorithm. 
Only the real part of the spectra is shown, and frequencies 
are relative to vh and vc- A positive or negative line in the 
spectrum indicates that the corresponding spin was in |0) or 
|1) before the read-out pulse. Furthermore, the position of 
the line in the spectrum of one spin also reveals the state of 
the other spin. For example, if the 1 H line is at vu — Jch/2, 
the 13 C spin is in |0); a 1 H line at vh + Jch/2 indicates the 
13 C spin is in Thus, the state of the two qubits for each 
of the four cases is (from top to bottom) |00), |01), 1 10) and 
jll). (Reproduced from Ref. [il] by permission of the Ameri- 
can Institute of Physics) 



FIG. 2. (a) Input and output states for the INEPT pulse 
sequence and (b) for the corresponding CNOT gate. 



FIG. 3. Schematic diagram of the main steps in quantum 
algorithms for period-finding. 



FIG. 4. (a) Spectrum of pentafluorobutadienyl cyclopen- 
tadienyldicarbonyliron complex in thermal equilibrium, (b) 
The same spectrum after preparing an effective pure state 
1 00000). (Reproduced from Ref. EH by permission of the 
American Physical Society) 
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